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ABSTRACT 

We present results based on Y JK^ photometry of star clusters in the Large Magel¬ 
lanic Cloud (LMC), distributed throughout the central part of the galaxy’s bar and 
the 30 Doradus region. We analysed the field-star decontaminated colour-magnitude 
diagrams of 313 clusters to estimate their reddening values and ages. The clusters are 
affected by a mean reddening of E{B — V)c. [0.2,0.3] mag, where the average internal 
LMC reddening amounts to ~ 0.1-0.2 mag. The region covering 30 Doradus includes 
clusters with reddening values in excess of E{B — V)= 0.4 mag. Our cluster sample 
spans the age range 7.0 ^ log(t yr“^) < 9.0, represents an increase of 30 per cent in 
terms of the number of clusters with robust age estimates and comprises a statistically 
complete sample in the LMC regions of interest here. The resulting cluster frequen¬ 
cies suggest that the outermost regions of the LMC bar first experienced enhanced 
cluster formation - log(t yr“^) G [8.5, 9.0] - before the activity proceeded, although in 
a patchy manner, to the innermost regions, for log(t yr“^) < 7.7. Cluster frequencies 
in the 30 Doradus region show that the area is dominated by very recent cluster for¬ 
mation. The derived star-formation frequencies suggest that the cluster and field-star 
populations do not seem to have fully evolved as fully coupled systems during the last 
- 100 Myr. 

Key words: techniques: photometric - galaxies: individual: LMC - Magellanic 
Clouds. 


1 INTRODUCTION 


* Based on observations obtained with VISTA at the Paranal 
Observatory under programme ID 179.B-2003. 
f E-mail: andres@oac.uncor.edu 


The galaxies in the Local Group are uniquely suitable to 
verify the results of analyses based on integrated cluster 


properties using resolved stellar photometry (e.g. de Grijs 


fc Anders|2006| [Colucci fc Bernstein|2012||Baumgardt et ah 


2013[ Cezario et al.||2013| de Meulenaer et al.||2013| ). Dedi- 

cated photometric surveys of the Magellanic Clouds at both 
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optical and near-infrared wavelengths are now probing suffi¬ 
ciently deeply so that we have a reasonable chance at re¬ 
solving individual stars to well below the main-sequence 
turn-off magnitudes characteristic of old stellar populations. 
With our unparallelled access to deep near-infrared obser¬ 
vations obtained with the 4 m Visible and Infrared Survey 


by the rapid removal of the intracluster gas on timescales of 


Telescope for Astronomy (VISTA; Emerson, McPherson fc 
Sutherland 2006[ Emerson & Sutherland 20101 as part of 


the VISTA near-infrared Y, J, Ka survey of the Magellanic 
System (VMC) the star cluster systems in the Small and 
Large Magellanic Clouds (SMC, LMC) can now be studied 
in detail to provide unique insights into the properties of 
their resolved star cluster populations. 

In this paper, we address a long-standing issue of con¬ 
tention in Magellanic Clouds studies, i.e., that of the cou¬ 
pling (if any) between field-star and (massive) star cluster 
formation (i.e., for cluster masses Md > 10'* M©). The LMC 
exhibits a well-known gap in the (massive) cluster age dis¬ 
tribution between ~ 3 and 13 Gyr, while the age distribu¬ 
tion of the field stellar population appears more continuous. 
Numerous authors have asserted that the LMC’s field-star 
and star cluster formation histories are significantly differ¬ 
ent (e.g. Olszewski, Suntzeff fc Mateo|1996 Geha et al. 1998 


Sarajedini 1998 and references therein). The situation for 
the SMC is less straightforward, although Rafelski & Zarit- 
sky (2005) provide tentative evidence that the cluster and 
field-star age distributions may also be significantly different 


in this system (see also Gieles, Lamers & Portegies Zwart 


2007). Similarly, the observed disparities between the clus¬ 
ter and field-star age distributions in the (dwarf) starburst 


galaxies NGC 1569 (Anders et al. 2004) and M82 (Barker, de 
Grijs fc Gervino|2008 ) seem to offer evidence in support of a 


decoupling between star cluster and field-star formation. If 
confirmed, this would be consistent with the view that mas¬ 
sive star cluster formation requires special conditions, e.g., 
large scale gas flows, in addition to the presence of dense 
gas (cf. Ashman fc Zepf||1992 Elmegreen fc Efremov|[l997 


de Grijs, O’Gonnell & Gallagher 2001). However, this sce¬ 


nario may well only apply to the most massive star clusters, 
whereas lower-mass clusters could indeed be following the 
field-star formation history more closely. This is our main 
scientific driver here. 

This is indeed an issue of general interest in the field 
of star and star cluster formation, given the commonly held 
notion that 70-90% of stars with masses in excess of 0.5 


Mq may form in clustered environments (e.g. Lada & Lada 


2003). Since most present-day stars in galaxies like the Milky 


Way and the Magellanic Clouds form part of their host 
galaxies’ field stellar population, it is not a stretch to suggest 
that a significant fraction of those field stellar populations 
originate from disrupted clusters, whether or not the latter 
were in fact gravitationally bound to begin with. 

There is now significant evidence that star cluster sys¬ 
tems appear to be affected by a disruption mechanism that 
acts on very short timescales (< 10-30 Myr) and which may 
be mass-independent, at least for masses > 10* Mq (e.g. 
Fall, Chandar fc Whitmore||2005l [Bastian et al.|[2005| |Fall| 
2006). This fast disruption mechanism, which is thought to 


disrupt up to 50-90% of the youngest clusters in a given 
cluster population (e.g. Lada fc Lada|199T Whitmore|2004 


Bastian et al.|20^5}|Mengel et al.|2005| [Goodwin fc Bastian| 

2006 Whitmore, Chandar & Fall|2007l) , is in essence caused 


< 30 Myr, a process coined cluster ‘infant mortality’ (Lada 


& Lada 2003); it was originally reported in the context of 


the number of very young embedded clusters in the Milky 
Way, compared to their older, largely gas-free counterparts. 

Those clusters that survive the infant mortality phase 
will be subject to the processes driving longer-term star 
cluster dissolution (see for reviews |de Grijs fc Parmentier 
2007 de Grijs||2010|. The longer-term dynamical evolution 


of star clusters is determined by a combination of internal 
and external timescales. These processes depend on the ac¬ 
tual cluster masses, with more massive clusters being more 
resilient to disruptive influences. Based on these arguments, 
one would expect that the lower-mass clusters in a given 
cluster population are the dominant donors of stars to their 
host galaxies’ field-star populations, while the higher-mass 
cluster formation and evolution history is not directly cou¬ 
pled to the field-star properties. In the context of the LMC’s 
star cluster population, cluster disruption has been the sub¬ 
ject of much debate. Nevertheless, and despite having access 
to cluster samples exceeding 900 members, the jury is still 
out as regards the typical time-scale of cluster disruption 


(e.g. Parmentier & de Grijs 2008) and whether this may vary 
as a function of position in the galaxy (e.g. [Bagheri, CionT 


|fc Napiwotzki|2013| |de Grijs, Goodwin fc Aiiders|2013[ and 

references therein). 

This paper is organised as follows. Section 2 provides 
an overview of the data on which this study is based and 
the data reduction procedures applied. We define our statis¬ 
tically complete LMG cluster sample in Section 3 and use 
isochrone fits to determine the clusters’ best-fitting phys¬ 
ical parameters in Section 4. Section 5 addresses the key 
questions posed in this study; in particular, we consider the 
cluster frequencies (GFs) across the LMG and compare the 
cluster formation rates with the star-formation rates (SFRs) 
in the corresponding field regions. We present our main con¬ 
clusions in Section 6. 


2 DATA HANDLING 

The VMC and its initial data are thoroughly described in 


Cioni et al. (2011), to which we refer the reader for details. 


Here we used data of two VMC tiles located in the central 
part of the LMC bar (LMC 6_4) and in a region covering 
the 30 Doradus star-forming region (LMC6_6), for which 
the VISTA imaging campaign has been completed. These re¬ 
gions are of significant scientific interest, since they contain 
a large number of star clusters which can provide statistical 
results about the central region of the LMC where other tiles 
are still being observed. Indeed, the purpose of this paper is 
to present a photometric analysis of the catalogued clusters 
located in those regions based on the most complete VMC 
data set to date. The available photometric data allow us to 
confirm the physical reality of the catalogued star clusters 
and estimate their fundamental parameters. We also investi¬ 
gate the CF in the regions of interest to assess whether star 
clusters and field stars have evolved as a coupled system. 
Figure 1 shows the location of all LMC tiles of the VMC 
survey. The red rectangles mark those used in this paper. 

The tile LMC 6_4 and 6_6 data refers to observations 
acquired from October 2010 until November 2011 and from 
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We performed PSF photometry using the 
iraf/daophot package to generate photometric cata¬ 
logues which return right ascension and declination, object 
magnitude, its error and sharpness, local magnitude com¬ 
pleteness and number of stars used to compute the latter in 
the Y, J, Ka passbands for ~2.6 million of sources in each 
tile. We used the psf task to produce the PSF model (which 
varies across the sky) and the ALLSTAR task to perform our 
photometry, using a radius of three pixels. We checked that 
our PSF photometry produced results consistent with those 
provided by the VSA for the bulk of the observed stars 
I Rubele et al. ]M2l |2015[ ). The final photometric catalogues 
are publicly available both at the VSAQ and ESC[^ science 
archives. 

We ran a large number of artificial-star tests to estimate 
the incompleteness and error distribution of our photometric 
catalogues for each tile and in every part of the colour- 


magnitude diagram (CMD), as described in Rubele et al. 
| |2012[ ) , where the photometry for two small subregions in tile 
LMC6_6 is presented. Photometric errors of 0.10 mag were 
derived for stars with Y = 20.8 mag, J = 20.6 mag and Ka 
= 19.7 mag in tile LMC 6_4 and for Y = 21.2 mag, J = 20.9 
mag, and Ka = 19.9 mag in tile LMC 6_6. The 50 per cent 
completeness level is reached at T = 19.9 mag, J = 19.5 mag 
and Ka = 19.6 mag in tile LMC 6_4 and at T = 20.4 mag, J 
= 20.1 mag, and Ka = 19.9 mag in tile LMC 6_6, throughout 
the entire tile areas. Photometric completeness is used here 
as a reference to our knowledge of the faintest cluster main- 
sequence turnoff (MSTO) reachable by our photometry. 


Figure 1. Spatial distribution of the |Bicaet al.| l |2008[ l’s catalogue 
of star clusters in the LMC centred at R.A. = 05“ 23™ 34°, Dec. = 
—69° 45' 22" (J2000), projected onto the sky. The objects studied 
in this paper are located in tiles highlighted as red rectangles: 
LMC6_4 (right) and LMC6_6 (left). The remaining VMC tiles 
across the LMC are shown as black rectangles. 


November 2009 to October 2010, respectively. Seeing con¬ 
straints, imposed for the purpose of homogenizing crowded 
and uncrowded field observations, range between 1.0 and 1.3 
arcsec (V), 0.9 and 1.2 arcsec (J), and 0.8 and 1.0 arcsec 
{Ka) and may exceed those values by 10 per cent accord¬ 
ing to observing policy. One tile covers uniformly an area 
of ~1.5 deg^, as a result of the mosaic of six paw-print im¬ 
ages, in a given waveband with 3 epochs at Y and J, and 12 
epochs at Ka. Individual epochs have exposure times of 800 
s (y and J) and 750 s {Ka). We used the vl.3 VMC release 
paw-prints. They were processed by the VISTA Data Flow 
System’s (VDFS; Emerson et al.|2004| pipeline (Irwin et al. 


2004[) and retrieved from the VISTA Science Archive (VSA; 


Hambly et al.|[2b04|. The processed paw-print images were 


used to derive the effective point spread functions (PSFs) 
using the iraf/daophot routines ( [Stetson, Davis fc Crab¬ 
tree ||1990|). We generated a reference PSF, which was con¬ 


volved with the paw-print images to homogenize the result¬ 
ing PSFs. We repeated these steps for each epoch separately. 
Finally, all homogenized paw-print images were combined 
using the SWARP tool ( Bertin et al.|2002 l, as described in 
[Rubele et al.| ( [2012[ ), thus generating deep tile images with 
homogeneous PSFs. 


3 THE CLUSTER SAMPLE 

Recognizing catalogued star clusters in deep VMC tile im¬ 
ages is neither straightforward nor simple. The catalogued 
objects were originally identified from optical images (e.g., 
from Digitized Sky Survey images; DS^^ which sometimes 
look rather different compared with their appearance at 
near-infrared wavelengths (bright red stars may dominate 
the cluster light while faint blue stars may be entirely ab¬ 
sent in the near-infrared). In addition, the spatial resolution 
and depth of the images on which the clusters were iden¬ 
tified differ from the equivalent parameters pertaining to 
the VMC images. Thus, for instance, single relatively bright 
stars might look like an unresolved compact cluster in im¬ 
ages of lower spatial resolution, or unresolved background 
galaxies could be mistaken for small star clusters in shal¬ 
lower images. Offsets in the compiled coordinates with re¬ 
spect to the objects centres cannot be ruled out either. 

To avoid mismatches between observed objects and the 
actual list of catalogued clusters, we first overplotted the 
positions of the catalogued clusters on the deepest stacked 


^ http://horus.roe.ac.uk/vsa/ 

^ http://www.eso.org/sci/observing/phase3/data_releases.html 
® The Digitized Sky Surveys were produced at the Space Tele¬ 
scope Science Institute under U.S. Government grant NAG W- 
2166. The images of these surveys are based on photographic data 
obtained using the Oschin Schmidt Telescope on Palomar Moun¬ 
tain and the UK Schmidt Telescope. The plates were processed 
into the present, compressed digital form with the permission of 
these institutions. 
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Ks image. This way, based on using the coordinates resolved 
by the SIMBAEj^astronomical data base, we unambiguously 
recognized the observed clusters one by one in the Ks image. 
Note that the main aim of this task is to confirm the com¬ 
piled cluster coordinates and sizes so as to extract from the 
VMC photometric catalogues the magnitudes of the stars in 
the cluster region. We are not interested in properties such 
as the clusters’ structures, stellar density prohles or radii, 
but we aim at obtaining accurate photometry for those stars 
which allow us to meaningfully dehne the clusters’ fiducial 
sequences in their CMDs. 

In order to build a cluster catalogue that is as complete 
as possible, we made use of those previously compiled by 


Pietrzyhski & Udalski (2000 hereafter PUOO), Bica et al. 
| 2008[ hereafer BOS) and Glatt, Grebel fc Koch| (2010 here¬ 
after GIO). Particularly, BOS provided angular dimensions, 
which we used for assessing whether such regions contain suf- 
hcient numbers of stars to construct reliable cluster CMDs 
with the smallest fraction of unavoidable field-star contami¬ 
nation. While recognizing the catalogued clusters on the Ks 
tile image, we discarded those in the line-of-sight of Hll re¬ 
gions, unresolved clusters, those for which the VMC photom¬ 
etry does not reach the MSTO or where the observed cluster 
CMDs do not show any cluster sequences. While attempting 
to cross-correlate the cluster names and their coordinates, 
we could not identify BSDL 1161 in the Ks tile image, and 
we found that SL36S is listed at incorrect coordinates in 
BOS (and, hence, in SIMBAD as well). Nevertheless, CIO 
compiled the correct cluster loci. Our final cluster sample 
includes 313 catalogued clusters. 


4 STAR CLUSTER PARAMETERS 

In general terms, the observed CMDs of the selected objects 
are the result of the superposition of different stellar pop¬ 
ulations distributed along the line of sight. For this reason, 
using the observed CMDs without subtracting the luminos¬ 
ity function and colour distribution of field stars may lead 
to incorrect interpretations. Moreover, since the catalogued 
clusters were initially identified as small concentrations of 
stars on the basis of stellar density fluctuations, their real 
physical nature requires subsequent confirmation. In order 
to disentangle cluster stars from field stars, we used CMDs 
of adjacent fields to subtract the local LMC field luminosity 
function and colour distribution. 

We constructed cluster CMDs based on all measured 
stars distributed within a circle with a radius of three times 
that tabulated by BOS based on visual inspection of the ob¬ 
jects on DSS images or by us from the deepest Ks images. 
The objects are of small angular size, typically ~ 0.6 arcmin 
(~ S.7 pc) in diameter. These regions are sufficiently large to 
encompass most of the cluster regions. Note that our main 
aim here is to clean the cluster CMDs from contamination 
by field stars within an area around the clusters’ centres that 
is nine times larger than Trr^iygjg,., so that we do not need 
to trace their radial profiles. Once the cluster areas were de¬ 
lineated, we defined four additional regions with areas equal 
to the cluster regions and located more or less equidistant 


* http://simbad.u-strasbg.fr/simbad/ 


from the clusters’ centres. For each of the latter regions we 
constructed CMDs showing the defining features of the local 
LMC star field in that particular direction. 

We next applied a procedure which was designed to 
compare each of the field CMDs to the cluster CMD and 
subtracted from the latter a representative field CMD in 
terms of stellar density, luminosity function and colour dis¬ 
tribution. This was done by comparing the numbers of stars 
counted in boxes distributed in a similar manner through¬ 
out all CMDs. The boxes were allowed to vary in size and 
position throughout the CMDs in order to meaningfully rep¬ 
resent the actual distribution of held stars. For additional 
details, see Piatti et al. (2014 2015aI, who applied this same 
procedure to other VMC data sets. Figure 2 illustrates the 
performance of the cleaning procedure for NCC 2108, where 
we plotted three different CMDs: that for the stars located 
within the cluster radius (top left-hand panel), a single-held 
CMD for an annulus - outer and inner radii equal to 3.163 
and 3.0 times the cluster radius - centred on the cluster, as 
well as the cleaned cluster CMD (bottom left-hand panel). 
The schematic diagram with a superimposed circle of radius 
equal to the cluster radius is shown in the bottom right- 
hand panel. The pink, light and dark blue hlled circles in 
the bottom panels represent stars with cluster membership 
probabilities P ^ 25 per cent, P = 50 per cent and P ^ 75 
per cent, respectively. The 313 individual photometric cat¬ 
alogues for the studied clusters are provided in the online 
version of the journal. The columns of each catalogue suc¬ 
cessively lists the magnitude, the error and the sharpness for 
every source within a circle of radius three times that of the 
cluster and centred on the cluster’s centre, in Y, J and Ks, 
respectively, the R.A. and Dec. and the photometric mem¬ 
bership probability (P). The latter is encoded with numbers 
1, 2, 3 and 4 to represent probabilities of 25, 50, 75 and 100 
per cent, respectively. 

We estimated reddening values and ages for our 313 
sample clusters using the theoretical isochrones of |Bressan| 
|et al.| ( |2012[ ) in the Vegamag system (where, by definition, 
Vega has a magnitude of zero in all filters), corrected by 
—0.074 mag in Y and —0.003 mag in Ks to put them on 
the VMC system ( [Rubele et al.|20T^ , to match the cleaned 
cluster CMDs. We adopted the same distance modulus for 


all clusters (m — M)o = 18.49 ± 0.09 mag (de Crijs, Wicker 
fc Bono|2014 1 and Ks — Mk^ = (m — M)q -|- 0.372P(P — V), 


for Rv = 3.1 ( Cardelli, Clayton fc Mathis] 1989| Cao 


et al. 20131, since by considering an average depth for the 
LMC disc of (3.44±1.16) kpc ([Subramanian fc Subrama- 


[200^ , we derived a smaller age difference than that 


resulting from the isochrones (characterized by the same 
metallicity) bracketing the observed cluster features in the 
CMD. 

We used isochrones for Z — 0.006 ([Fe/H] = —0.4 dex), 
which corresponds to the mean LMC cluster metal content 
for the last ~ 2-3 Cyr ( Piatti fc Ceisler|2013[ ). Additionally, 
the Y — Ks colour is not sensitive to metallicity differences 
smaller than A[Fe/H] ~ 0.4 dex, which is adequate given the 
spread of the stars in the CMDs (see Piatti et al.||2015a| ). 
Indeed, we tried using isochrones with metallicities [Fe/H] 
= 0.0 dex and [Fe/H] = —0.7 dex and found negligible dif¬ 
ferences with respect to that of [Fe/H] = —0.4 dex, keeping 
in mind the relatively sparse nature of the majority of our 
clusters and the intrinsic spread of the stars in the CMDs. 
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Table 1. Fundamental parameters of the LMC clusters studied 
here. We assume a distance modulus of 18.49 mag and a metal- 
licity of [Fe/H]= —0.4 dex. The complete table is available online 
as Supporting Information. 



Figure 2. CMDs for stars in the field of NGC 2108: the observed 
CMD composed of the stars distributed within the cluster radius 
(top left-hand panel); a field CMD for an annulus - outer and in¬ 
ner radii equal to 3.163 and 3.0 times the cluster radius — centred 
on the cluster (top right-hand panel); the cleaned cluster CMD 
(bottom left). Colour-scaled symbols represent stars that statisti¬ 
cally belong to the field (P Si 25%, pink), stars that might belong 
to either the field or the cluster (P = 50%, light blue), and stars 
that predominantly populat e the cluster region ( P ^ 75%, dark 
blue). Three isochrones from |Bressan et al.| | 201^ for log(t yr“^) 
= 8.8, 8.9, and 9.0 and Z = 0.006 are also superimposed. The lack 
of fainter main-sequence stars results in less satisfactory match¬ 
ing at fainter magnitudes. The schematic diagram centred on the 
cluster for a circle of radius three times the cluster radius is shown 
in the bottom right-hand panel. The black circle represents the 
adopted cluster radius. Symbols are as in the bottom left-hand 
panel, with sizes proportional to the stellar brightnesses. North 
is up; East is to the left. 


We found that isochrones bracketing the derived mean age 
by AE{B — V) = ±0.05 mag and A log(t yr“^) = ±0.1 rep¬ 
resent the overall uncertainties associated with the observed 
dispersion in the cluster CMDs. Although the dispersion is 
smaller in some cases, we prefer to retain the former val¬ 
ues as an upper limit to our error budget. Table 1 lists the 
diameters used and the resulting reddening and age values 
obtained for the entire cluster sample, while the bottom left- 
hand panel of Fig. 2 illustrates the results of the isochrone 
matching for NGC 2108. 

The resulting reddening distributions for the clusters 
located in tiles LMC 6_4 and 6_6 are shown in Fig. 3, where 


Name 

R.A. 

(°) 

Dec. 

n 

d 

(arcmin) 

E{B - V) 
(mag) 

log(t yr 1) 

NCC 1858 

77.500 

-68.904 

0.70 

0.30 

7.30 

HS335 

82.737 

-69.341 

0.40 

0.10 

8.60 

NGC 2009 

82.750 

-69.182 

0.80 

0.15 

7.50 


we have also included an inset showing their spatial dis¬ 
tribution. We converted E{Y — Kb) to E{B — V) using 
E{Y - Kb) = 0.84:E{B - V) ( jCardelli, Clayton fc Mathis 
19891. Both surveyed regions are characterized by a mean 


reddening spanning the range E{B — V) £ [0.2,0.3] mag. 
The tile centred on the LMC bar contains a slightly larger 
number of clusters with reddening values towards the up¬ 
per limit of this range. On the other hand, the tile cover¬ 
ing 30 Doradus includes clusters with reddening values in 
excess of E{B — V) = 0.4 mag. Bearing in mind that the 
observed cluster sequences are affected by interstellar extinc¬ 
tion which is the sum of Milky Way (foreground) and LMC 
(internal) reddening, and that the mean Galactic reddening 
towards the LMC is E{B — F) « 0.08 mag (Dolphin et al. 
2001[ Schlafly fc Finkbeiner|[2bll I, it follows that the aver¬ 
age (internal) reddening in these particular LMC regions is 
0.1-0.2 mag. This range is in very good agreement with 


a number of previous estimates (e.g.. 

Oestreicher, Cocher- 

mann & Schmidt-Kaler 

1995||Subramanian Subramaniam 

2009| Tatton et al.|2013 

1, although slightly lower values have 


also been found from analysis of field red-clump and varl 
able stars (Clementini et al.|2003 Haschke, Grebel & Duffau 


20111 . 


Figure 4 shows a comparison of the cluster age estimates 
derived here with those previously obtained from similar 
isochrone matching. The latter assumed a constant metal 
content for the LMC of [Fe/H] = —0.4 dex, a distance mod¬ 
ulus of 18.50 mag (except PUOO who adopted 18.24 mag and 
mentioned that change in the distance modulus by about 
0.15 mag produced difference in ages of less than a few 
percent, which is much less than the derived errors), and 
isochrones computed by the Padova group ( [Bertelli et al.j 
1994 Marigo et al.|20d8 |. Agreement with the OGLE-based 
age estimates (200 clusters) is good, just like that with the 
cluster ages determined on the basis of Washington pho¬ 
tometry (Piatti 2012[ 2014a [Choudhury, Subramaniam fc 
Piatti 2015 3, 7 and 3 clusters, respectively). The most 


discrepant points in the OGLE age comparison are associ¬ 
ated with clusters whose CMDs include relatively bright blue 
stars that PUOO either rejected or included in the isochrone 
matching, thus resulting in, respectively, older (HS 223, log(t 
yr“^) = 8.9; BSDL 725, log(t yr“^) = 9.0) or younger 
(HS362, log(t yr-i) = 6.9; HS219, log(t yr"^) = 7.0; H88- 
169, log(t yr“^) = 7.1) ages than our values. As for the 
comparison with the ages based on Washington photome¬ 
try, which were obtained from isochrone matching to CMDs 
previously corrected for field-star contamination using the 
same technique described above and for clusters not stud- 
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Figure 3. Reddening histograms for tiles LMC6_4 (black) and LMC6_6 (red) (bottom panel) and the spatial distribution of the clusters’ 
reddening, with darker shading representing higher reddening values (top panel). The grey scale is in units of E{B — V) (mag). 
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Figure 5. Spatial distribution of the cluster ages, with darker shading representing older ages, centred at R.A. = 05^ 23*^ 34®, Dec. = 
—69° 45^ 22^^ (J2000). North is up; East is to the left. The grey scale is in units of log(t yr“^). The locations of 30 Doradus and the 
LMC bar are indicated with a dashed circle and two parallel lines ( [Subramaniam Sz Subramanian|2009[ |, respectively. 
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Figure 4. Comparison of the cluster age estimates derived here 
with those from the literature: PUOO (top panel), GIO (middle 
panel) and Washington photometry | |Piatti|2012| |2014a[ |Choud-| 
|hury, Subramaniam fc Piatti||2015| l (bottom panel). Error bars 
and the locus of equality are also shown. 


led by either PUOO or GIO, better agreement can still be 
achieved. The correlation with the ages determined by GIO 
(89 clusters) is less tight, because GIO did not perform field- 
star decontamination, although these authors mention that 
held contamination is a severe effect in the extracted cluster 
GMDs and may, therefore, affect the age estimates signif¬ 
icantly. Consequently, their age estimates likely rehect the 
composite LMC stellar populations (e.g., KMK 88-74 (log(t 
yr"^) = 7.8), BSDL2807, (log(t yr-^) = 8.0), KMK88-86 
(log(f yr”'^) = 8.4)). This possibility alerts us to the fact 
that solely the circular extraction of the observed GMDs of 
clusters located in highly populated stellar helds is not suf- 
hcient for accurate isochrone matching of the main cluster 


features (see also 

Piatti & Bica 2012 

Piatti 2014a Piatti 

et al.|2014| 2015a 

I. 


Based on Fig. 4, we conclude that the age scale de¬ 
rived in this paper is consistent with those pertaining to our 
comparison samples. We can therefore safely merge these 
samples to compile a larger cluster sample characterized by 
an age distribution put on an homogeneous scale. Of our 
313 sample clusters, 95 did not have any previous age es¬ 
timates - a ~ 30 per cent increase of clusters with known 
ages in these fields. The final cluster age catalogue is com¬ 
posed of our 313 age estimates, combined with 74 clusters 


from PUOO, five clusters from GIO and 15 clusters from other 
studies (mainly based on Washington photometry). When¬ 
ever a cluster has more than one age estimate, we adopted a 
weighed average value. All ages considered here were derived 
based on isochrone matching of the cluster GMDs, spanning 
the age range 7.0 < (log(t yr“^) < 9.0. Figure 5 shows the 
clusters’ spatial distribution, with darker points represent¬ 
ing older ages. Note the obvious distribution of very young 
clusters around 30 Doradus in tile LMC 6_6, as well as of 
young clusters in the innermost regions of the bar in tile 
LMC 6_4. 


5 THE CLUSTER FREQUENCIES IN TILES 
LMC 6_4 AND 6_6 


A number of different studies have shown that the LMC’s 
held SFR is variable across the galaxy ([Rubele et al.||2012 


[Meschin et al.| |2014[ [Harris fc Zaritsky| |2009| hereafter 
HZ09). In particular, HZ09 concluded from the concordance 
between the star-formation and chemical-enrichment histo¬ 
ries of the held and cluster populations that the held and 
cluster star-formation modes are tightly coupled. Recently, 


Piatti & Geisler (20131 studied 21 LMC regions spread 


across the main body of the galaxy and found that the clus¬ 
ter and held star age-metallicity relationships (AMRs) of 
such regions show a satisfactory match only for the last 3 
Gyr, i.e., log(t yr“^) ^ 9.5, while for the oldest ages - t >11 
Gyr, log(t yr“^) > 10.05 - the cluster AMR represents a 
lower envelope to the held AMR. 

We are interested in constructing the CFs, i.e., the num¬ 
ber of clusters formed per unit time interval as a function of 
age, as a function of position in both tiles studied here. It 
would appear reasonable to infer that if the star-formation 
history of LMC held stars has been different throughout the 
galaxy, and held stars and clusters show some evidence of 
common formation and chemical-enrichment histories, then 
the CFs should rehect the same spatial variations as seen 


in the held stars. Piatti (2014b hereafter P14) showed that 
the number of catalogued clusters without age estimates is 
relatively small with respect to the total number of cata¬ 
logued clusters in the LMC. He also showed that nearly half 
of the more massive clusters in his sample with age estimates 
mostly trace the overall behaviour of the CFs. Alternatively, 
clusters with masses below 10^ Mq for ages younger than 1 
Gyr do not contribute signihcantly to the normalized CFs, 
although they represent on average nearly half of the total 
number of clusters with age estimates. In addition, he found 
that different model age distributions for the total number 
of clusters without age estimates do not affect the resulting 
CFs. Moreover, he concluded that possible incompleteness 
of non-catalogued clusters (which include the faintest clus¬ 
ters in the galaxy) does not play an important role in the 
resulting CFs either. Since our cluster age catalogue for the 
two LMC tiles of interest contains ~ 30 per cent more clus¬ 
ters than previously known, we are confident that we can 
consider our results statistically robust. Indeed, our cluster 
sample includes ~ 80 per cent of the catalogued clusters in 
the surveyed regions. 

In constructing the CFs, we first divided the tiles into 
12 subregions of similar size, and then counted the number 
of clusters in each region, taking into account bin size and 
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age uncertainty effects. This procedure allowed us to build 
age histograms (and hence CFs) which best reproduce the 
intrinsic age distribution of each subregion. For additional 
details about the treatment of bin size and age uncertainty 


effects, see, e.g., Piatti & Geisler (20131, P14 and Piatti 


et al. |2015a|. All CFs were normalized to the total number 


of clusters included in our sample. We assume that there 
has been no statistically significant mixing of clusters for 
ages of log(t yr“^) < 9.0) (which represents less than the 
last 10 per cent of the galaxy’s lifetime) between regions. 
Indeed, |van der Marel fc Kallivayalil (20141 showed that for 
most tracers (including clusters and field stars), the LMC’s 
line-of-sight velocity dispersion is at least a factor of ~2 
smaller than the corresponding rotation velocity, which im¬ 
plies that the LMC as a whole is a kinematically cold disc 
system. Figures 6 and 7 show the resulting CFs for the sub- 
regions of tiles LMC 6_4 and 6 _ 6 , respectively. We include 
insets with the schematic distribution of the 12 subregions 
and a colour reference of the respective CF curves. From Fig. 
6 , it is clear that the outermost subregions of the LMC bar, 
which crosses tile LMC 6_4 along the SE-NW diagonal (blue 
lines), have experienced more prominent cluster formation 
in the past - log(t yr“^) ~ 8 .5-9.0 - and more recently, log(t 
yr“^) < 7.7, cluster formation in the innermost regions has 
become more active (red and green lines). However, the clus¬ 
ter formation scenario in the central parts of the bar appears 
to be rather complex and patchy. In tile LMC 6_6 (Fig. 7) 
the most striking role is played by the 30 Doradus region, 
which dominates the recent cluster formation activity (dark 
green and black lines). Likewise, 30 Doradus is spatially sur¬ 
rounded by cluster formation activity characterized by older 
ages with increasing distance from the object (blue and red 
lines), possibly approaching the intrinsic CF in the outer 
LMC bar (see fig. 5 in P14). 


We also compiled composite CFs using all clusters lo¬ 
cated within each tile and compared them with those previ¬ 
ously obtained by P14 for HZ09’s bar and the 30 Doradus 
region (see hg. 6 in HZ09). The CFs were constructed by 
taking into account bin sizes and age uncertainties, and nor¬ 
malized to the total number of clusters. Figures 8 and 9 
show the resulting CFs and those from P14, represented by 
black and grey lines, respectively. Both pairs of CF curves 
are very similar. CF differences become important for log(t 
yr“^) G [7.2, 7.6]. These differences are likely owing to the 
fact that P14’s CF covers the entire LMC bar, while the 
present CF targets its central part only, which is charac¬ 
terized by a relative excess of younger clusters with respect 
to the entire bar. On the other hand, the most noticeable 
differences between the CFs centred on 30 Doradus occur 
for older ages - log(t yr“^) > 8.5 - which is most likely 
also caused by the larger area covered by HZ09 than our 
tile LMC 6 _ 6 , which thus implies that the HZ09 CF is com¬ 
posed of a relatively smaller fraction of older clusters. We 
superimposed on Figs. 8 and 9 the CFs (red lines) and the 
corresponding uncertainties for the bar region in HZ09 as 
well as their 30 Doradus region, respectively, assuming that 
clusters are born according to a power-law mass distribution 
with a slope a = —2 and at a rate that is proportional to the 
SFR determined by HZ09 for each region. The shapes of the 
CFs recovered from the SFR curves generally follow those of 
the observed CFs, although the cluster and field-star popu¬ 
lations do not seem to have evolved as fully coupled systems. 


Nevertheless, coupling between field-star and cluster forma¬ 
tion appears to have been dominant for both regions. 

Assuming that most - if not all - stars form in clus¬ 
tered modes (Lada & Lada 20031, our results suggest that 
the combination of cluster formation and disruption rates 
(which depend on environmental conditions, cluster masses, 
among others) still provides a positive balance, in the sense 
that the cluster population has not been entirely trans¬ 
formed into field stars. [Baumgardt et al.| ( [2013| ) showed that 
about 90% of all clusters older than 200 Myr are lost per 
dex of lifetime, which implies a cluster dissolution rate sig- 
nihcantly faster than that based on analytic estimates and 
N-body simulations. However, |de Crijs, Coodwin fc Anders| 
(20131 showed that there is no evidence of significant de¬ 
struction, other than that expected from stellar dynamics 
and evolution in simple population models for ages up to 1 
Cyr (log(t yr“^) = 9). As can be seen, the issue about at 
what pace cluster disruption takes place in different LMC 
regions and even whether field stars mostly come from clus¬ 
ter disruption still needs further study. Our results do not 
allow us to draw any further conclusion about this issue, but 
that the cluster formation has been different in the studied 
regions. 

By comparing our results with those of HZ09 for the 
two tiles, some hint about the distint roles played by the en¬ 
vironments can be drawn, since differences in CF are larger 
in the tile LMC 6 -6 than in LMC 6-4 for ages younger than 
log(t yr“^) < 8.0. Werchan & Zaritsky (2011 1 tabulated the 
concentration, central surface brightness, tidal radii, 90% 
enclosed luminosity radii, and local background luminos¬ 
ity density for 1066 star clusters in the Magellanic Cloud 
Phomotometric Survey (MCPS [Zaritsky et al.|[2002[ ). They 
found that while most of the star clusters are similar be¬ 
tween the LMC and the SMC populations, the LMC lacks 
star clusters that are large, either in terms of core or 90% 
enclosed luminosity radii as the largest in the SMC, and 
suggested that this result could be featuring a signature of 
larger tidal stresses in the LMC. On the other hand, [Wilkin- 1 


son et al. (20031 used N-body simulations of star clusters to 


investigate the possible dynamical origins of the observed 
spread in core radius among intermediate-age and old star 
clusters in the LMC and concluded that the tidal field of 
the LMC has not yet signihcantly influenced the evolution 
of the intermediate-age clusters. In this sense, we interpret 
that the actual dynamical state of the studied LMC clusters 
might be a combination of an age effect (two-body relatax- 
ion) and their locations in the galaxy (see., e.g. [Clatt et al. 


2011 for the SMC). 


Finally, we compared the CFs derived here with those 
obtained from the SFRs derived by Rubele et al. (20121 for 
tile LMC 6 - 6 . The latter authors derived the SFRs for two 
subregions (D1 and D2; see their fig. 2). They used the ‘par¬ 
tial model’ technique, which applies a linear combination 
of partial models that optimally matches the observed Hess 
diagrams allowing the distance modulus, visual extinction, 
age and metallicty to vary, convolved with the distributions 
of photometric errors and completeness. The coefhcients of 
this linear combination of partial models (including the best- 
htting distance modulus and extinction) are directly con¬ 
verted into the SFRs. We used their SFRs and cluster masses 
from log(Mci[M 0 ]) = 2.2 to log(Mci[M 0 ]) = 5.0 ( jde Crijs~fc| 
Coodwin|2008 [Clatt et al.|20li|, and normalized the result- 
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log(t) [yr-i] 


Figure 6. CFs for subregions in tile LMC 6_4. The bottom right-hand panel shows the spatial distribution of |Bica et al. ]l |2008l l ’s catalogue 
of star clusters, while the left-hand panel schematically represents the subregions and indicates the colour used in drawing the respective 
CFs. North is up and east is to the left. Note that computed values (filled circles) for each CF are connected by a solid line. Some CFs 
do not have points for certain age intervals because of the lack of clusters within them, which resulted in truncated CFs. In such cases, 
we drew isolated points with a short line to represent that they belong to a truncated CF. One subregion does not contain any cluster 
in our sample. 
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Figure 7. As Fig. 6, but for tile LMC6_6. 


ing CFs by the total number of clusters used, so that they 
can be compared directly with the observed distributions. 
Figure 10 shows the resulting, recovered CFs, where the un¬ 
certainties are indicated by dashed lines. The observed and 
recovered CFs are clearly different for ages younger than 


log(t yr“^) ~ 8.0. Even though the recovered CFs require 
additional refinements, the cluster excesses could be evi¬ 
dence of relatively recent enhanced cluster formation in this 
part of the galaxy. Likewise, note that both recovered CFs 
are different between them, which suggests either that clus- 
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log(t) [yr-i] log(t) [yr-i] 


Figure 8. Composite CF for tile LMC6_4 (black) compared with 
those obtained by P14 for clusters within HZ09’s bar region (grey) 
and by HZ09 for the respective mean SFR and its Itr uncertainty 
(solid and dashed red lines, respectively). 


ter disruption changed dramatically from one subregion to 
the other, or field stars did not come from cluster disrup¬ 


tion in a similar proportion. Rubele et al. (2012 \ and Rubele 


et al. (20151 showed that the SFR of field stars varies with 


the position in the LMC and SMC, respectively. Conversely, 
P14 and Piatti et al. 12015b I found evidence that a similar 
behaviour occurs within the LMC and SMC cluster popu¬ 
lations, but not with a completely equivalent pace to field 
stars. These results, besides those derived here, provide some 
clues for a better approach in the study of the field stars 
origin and its link to cluster disruption and environmental 
conditions. 


6 CONCLUSIONS 

In this paper, we have analysed CMDs of previously cata¬ 
logued star clusters located in the innermost region of the 
LMC based on a YJKa photometric data set obtained by 
the VMC collaboration. We focused on tiles LMC 6_4 (lo¬ 
cated in the central part of the LMC bar) and 6_6 (centred 
on 30 Doradus), because they are among the first fully com¬ 
pleted tiles of the VMC survey, and they show the most 
recent star-formation activity in the galaxy. We obtained 
PSF photometry for stars in and projected along the line 


Figure 9. As Fig. 8, but for tile LMC 6_6 and HZ09’s 30 Doradus 
region. 


of sight towards 313 catalogued clusters, which represents a 
statistically meaningful sample size. 

We applied a field-star subtraction procedure to statis¬ 
tically clean the cluster CMDs from field-star contamination 
and to disentangle cluster features from those associated 
with their surrounding fields. The technique we employed 
makes use of variable cells to reproduce the field CMDs as 
closely as possible. Based on matching theoretical isochrones 
in the VISTA system to the cleaned cluster CMDs, we ob¬ 
tained estimates of the clusters’ reddening values and ages, 
assuming a distance modulus of 18.49 mag and a metallic- 
ity of Z = 0.006 ([Fe/H] = —0.4 dex). Both surveyed re¬ 
gions are affected by a mean reddening spanning the range 
E{B — V)€ [0.2, 0.3] mag; the average internal (LMC) red¬ 
dening amounts to ~ 0.1-0.2 mag. In addition, the tile cov¬ 
ering 30 Doradus includes clusters with reddening values in 
excess of E{B — V) — 0.4 mag. 

The resulting cluster ages span the age range 7.0 ^ log(t 
yr“^) < 9.0. They form part of the cluster data base which 
will result from the VMC survey and which will be used to 
self-consistently study the overall cluster-formation history 
of the Magellanic system. The age scale derived in this paper 
is consistent with those resulting from previous publications. 


in particular those based on OGLE (Pietrzyhski & Udalski 


20001, the Magellanic Cloud Photometric Surveys (Zaritsky 
et al.||2002[) and the Washington photometric system (|P1 


atti|2012 2014a Choudhury, Subramaniam fc Piatti|20i5 1, 
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log(t) [yr-i] 


Figure 10. CFs for subregions D1 (thick red line) and D2 (thick 
blue line) in tile LMC 6_6. Thin and dashed lines represent, re¬ 
spectively, the mean CFs and their Itr dispersions recovered from 
the corresponding SFRs obtained by |Rubele et ah'] | |2012[ |. 

among others. A total of 95 clusters of the 313 objects in our 
sample had their ages determined for the first time. This rep¬ 
resents an increase of ~ 30 per cent of clusters with known 
ages in the two LMC tile regions. The enlarged catalogue of 
clusters with age estimates contains 407 objects. 

Based on analysis of this statistically complete cluster 
sample, we addressed the nature of variations in the LMC 
CF in terms of the clusters’ spatial distribution. Particu¬ 
lar attention was paid to bin sizes and age uncertainties 
in constructing CFs for 12 different subregions of similar 
sizes, homogeneously distributed within each LMC tile. As 
for tile LMC 6_4, we found that the outermost subregions 
of the LMC bar experienced enhanced cluster formation ac¬ 
tivity in the past - log(t yr”'^) G [8.5,9.0] - whereas more 
recently, log(t yr“^) < 7.7, cluster formation seems to occur 
more frequently in the inner regions. Overall, however, clus¬ 
ter formation in the central parts of the bar has proceeded 
in a patchy manner. On the other hand, in tile LMC 6_6, 
the 30 Doradus subregions play the most striking role in the 
very recent cluster formation history. They are surrounded 
by successively older clusters for increasing distances from 
30 Doradus. 

Finally, we compared the composite CFs for each LMC 
tile to those recovered assuming that clusters have formation 
rates similar to the known field SFRs. The shapes of the re¬ 


covered CFs, and more particularly those for the subregions 
near 30 Doradus, suggest that cluster and field-star popula¬ 
tions do not seem to have evolved as fully coupled systems 
during the last ~ 100 Myr. The cluster excesses found could 
be evidence of relatively recent enhanced cluster formation 
in these parts of the galaxy and thus confirm both that most 
stars are formed in clusters and that there is no significant 
destruction other than that expected from stellar dynamics 
and evolution of simple stellar population models for ages 
up to 1 Gyr. 
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